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We consider a model Hamiltonian for a dimer of length a including all the electronic one- and two- 
body terms consistent with a single orbital per site, a free Einstein phonon term for a frequency fi, 
and an electron-phonon coupling go of the Holstein type. The bare electronic interaction parameters 
were evaluated in terms of Wannier functions built from Gaussian atomic orbitals. An effective 
polaronic Ifamiltonian was obtained by an unrestricted displaced-oscillator transformation, followed 
^ , by evaluation of the phononic terms over a squeezed-phonon variational wave function. For the cases 

of quarter-filled and half-filled orbital, and over a range of dimer length values, the ground state for 
given (?o and Q. was identified by simultaneously and independently optimizing the orbital shape, the 
^ , phonon displacement and the squeezing effect strength. As a varies, we generally find discontinuous 

changes of both electronic and phononic states, accompanied by an appreciable renormalization of 
the effective electronic interactions across the transitions, due to the equilibrium shape of the wave 
functions strongly depending on the phononic regime and on the type of ground state. 
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I. INTRODUCTION 



> . 

■^■^ . Much attention has been devoted to materialspwhere structural effects deep y influence the electronic phenom- 
00 ' ena, as, for instance, V2-a;Cra:03,EJ polyacetilenep^El displacive-type ferroele(fijrics,cl the manganites exhibiting colossal 
r magnetoresistancej3 and the high temperature superconducting perovskites.QTl From the theoretical standpoint, such 
' effects have been studied by adding various kinds of electron-phonon couplings to different electronic Hamiltonians, 
\ mainly of the extended Hubbard or t-J type. However, comparatively little attention has been paid to the effect of 
the structure (in its simplest form, a variation of the lattice spacing a) onto the electronic interactions themselves. 
Such effects are of great relevance in experiments: consider, for example, the measurements under cither external 
I or chemical pressure, and the wealth of interesting phenomena that they have revealed. The electronic interaction 
parameters can be written as integrals of various combinations of Wannier functions centered on the lattice sites (as 
I will be detailed in Section |l| below). Therefore, a change of a is expected to have a non-negligibJe effect on them. 
'■^ In the pe^lm of electronic transitions, to our knowledge, after a pioneering attempt by KawabataEl, only Spalek and 
^ , W6jcik£2l focused onto such effects. They argued that, in a system of fermions undergoing a metal-insulator transition 
O ■ (MIT) due to the on-site repulsion effect, the wave functions may significantly change their shape across the transition, 
concomitantly to a change of the equilibrium lattice spacing. Therefore, different values for the hopping amplitude 
tij and the on-site repulsion U may correspond to the itinerant and the localized regimes. In the present work we 
explore the consequences of such an approach in systems with a Holstein-type electron-phonon interaction. Our model 
Hamiltonian includes all the electronic interaction parameters for one- and two-body terms, whose dependence on 
a results from their evaluation in terms of Wannier functions built from Gaussian orbitals. We also consider the 
characteristic structural length a as free to vary between reasonable boundaries, in order to model the effects of the 
experimental techniques that allow for a more or less continuous tuning of the bond length. We find that both the 
electronic and phononic states significantly depend on a. For instance, on one side, the phonon state, characterized 
by the extent of the displacement and the squeezing, reflects the electronic transitions at a deeper level than just a 
frequency renormalization. On the other side, as the optimal shape of the Wannier wave functions depends not only 
on a, but also on the phonon state, the renormalization of the electronie-ojarameters by the phonons is actually more 
complex than indicated by the standard results of the polaronic model.til 

To put into evidence the basic features of our approach, avoiding unnecessary complications, our model system in 
the present paper is a simple dimer with a single orbital per site. It can be considered as the basic unit to build 
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clusters and lattices, and its simplicity allows for a full control on the analytic aspects of the calculations. Also, 
some unavoidable approximations which have to be introduced are linked to the physics, and not to the complications 
of treating larger systems. Some new aspects distinguish the present work from the rich existing literature on the 
Holstein dimer with various types of inter-electronic interactions, (see, for instance, Refs. 3, |l2| , ^ and the literature 
cited therein). First, the effective electronic Hamiltonian originates from a first- principles Hamiltonian including the 
complete set of one- and two-body electronic terms consistent with a single orbital per site. To the best of our 
knowledge, previous studies considered only a more restricted subset of interactions. Second, we present our results 
mainly by setting the values of the electron-phonon coupling go and phonon energy hfl and taking the dimer length 
a as the continuously varying controlling parameter. The dependence of the state of the system on its length usually 
receives little attention even though it is much more easily accessible experimentally than the dependence on the 
strength of the coupling, which requires delicate chemical manipulations. However, some discussion of the case of 
fixed length and varying coupling will also be given. We will discuss the cases of quarter- and half-filling of the electron 
orbitals. 

The paper is organized as follows: Section || presents the model Hamiltonian and the expressions of the electronic 
interactions explicitly depending on the dimer length. In Section [H the effective electronic Hamiltonian is obtained, 
while in Section IV several correlation functions and other quantities of experimental interest are explicitly evaluated. 
The results of the numerical analysis of our model are presented and commented in Section 0. Finally, some general 
conclusions are presented in Section VI. 



II. THE MODEL 

For a dimer consisting of identical ions centered at Ri and R2 we consider the Hamiltonian H = H^^ + H^^ + H^^ 
where Hd includes all the electronic one- and two-body terms for a single orbital, H^^ is the free-phonon term for a 
single Einstein frequency fi, and H^p introduces an electron-phonon coupling qq of the Holstein type. 

In standard notation, we have 



Vn,n, - 2J,SIS'^ - J.,y{S+S- + H.c.) + P{c\^c\^c^^c^^ + H 
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(1) 
using 



The bare electronic parameters eo, t, X, U, V, and P {= Jz — Jxy) were evaluated according to Ref. 
normalized Gaussian "orbitals" 0i(r) = (2r^/7r)'^^^ exp [— r^(r — R^)^] {i = 1,2). Their overlap integral is S' = 
{4'i\4'2) = exp(— r^a^/2), where a — |Ri — R2I is the length of the dimer. Then the two orthonormalized Wannier 
functions are 



*i(r) = A+(5)0i(r)+A_(5)02(r) 
*2(r) =A_(5)0i(r)+A+(5)02(r) 



with 



1 



± 



1 



(2) 



(3) 



obtain 



r; I ) [1 + Fo(2a2r2) - 2S'Fo{a'T' /2)\ , 



In order to calculate the parameters one has to evaluate integrals involving Gaussian functions, which have been 
discussed in Ref. ^ If the ion charge is Z\e\, and defining /3 = 7i^/2m, rj = — Ze^, and Fo{x) = x~^/^Erf (x^^^), we 

(4) 

(5) 

(6) 
(7) 




[1 + 252 + Fo(r'a') - 2(1 + S^)Fo{r^ay4)] , 
[2-S^ + 25"* + S^Fo{T^a^) - AS^FoiT^a^ /A)] , 
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(1-^2)2 



[5^(1 + 25-2) + (2 - S^)Fq{T'^c?) - 4S'2Fo(r2aV4)] , 



(1 -S'2)2 



[3 + Fo(r2a2)-4^^o(r2a2/4)] 



(8) 
(9) 



The parts of the Hamihonian involving phonons are a free-phonon term and a Holstein-type interaction iJop, 
reading 



and 



(10) 



(11) 



where M is the mass of the ions, is the phonon frequency, assumed dispersionless for simphcity (Einstein phonon), 
and Pi = —ihd/dui is the momentum of the ion at position i. We are considering only longitudinal vibrations of the 
dimer. 

The standard procedure to quantize the phononic Hamiltonian requires the introduction of usual phonon operators 

a^ = {2Mhny^/^ {Mflu, + iP^) (12) 
for each site. Actually, it is more convenient to work with symmetric (s) and antisymmetric (&) combinations 

s = (fli + a2)/\/2 , 



6 = (fli - a2)/\/2 , 
such that the free-phonon Hamiltonian H^^^ can be decomposed as 

^ ^ph + ^p°h = (b'^b +^^+hn[ 

and the electron-phonon Hamiltonian iJop has two contributions 

H. 



- hnjo{ni + ^2) (s^^ + s) , 
hn-foim - n2) {b'' + b) , 



(13) 



(14) 



(15) 



where we have used the notation rii = X^o-'^io-, and introduced the coupling energy per unit charge ga = g^^/%/2MQ, 
which, conveniently scaled by the phonon energy, yields the dimensionless quantity 70 = go/hfl. Eqs. (15) show 
explicitly the coupling of s with the total charge, and b with the charge transfer. 



III. THE EFFECTIVE ELECTRONIC HAMILTONIAN 

Following the procedure of Ref. ^ we will obtain the effective electronic (polaronic) Hamiltonian Hpo\ by first 
performing on H a, generalized Lang-Firsov (displaced-oscillator) transformation. Then, in the transformed Hamilto- 
nian, we will take the average of the phononic parts over the squeezed-phonon wave function to eliminate the phonon 
degrees of freedom. 

As we have two different boson operators, s and b, there will be two corresponding generators Rs and i?f,, and we 
have to perform a sequence of two unitary transformations on H. The generators are defined as 



^^6 = ^7o("i - {b^ - b) 



(16) 



and 



Rs ^ Vlnini + n2) [s"^ - s) , (17) 
where 5 and 77 are so far undetermined parameters, whose value will be fixed later in a convenient way. 
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Notice that [Rb, Rs] = 0. Then, the transformed phonon operators are 

e^'be-"" =b- 5-/o{ni - ^2) (18) 



and 



We transform, therefore, H^^ + H^^ as 



e =se ^ = s — 7y7o(ni + . (19) 



(^b^b + + - S)jo{ni - n2){b^ + b) 

nn5{2 - 5)-il{ni - n2f (20) 



and 



2, 

- nnTji2 ~ Tj)j^{ni + n2 f . (21) 

Since the purely electronic terms of the total Hamiltonian commute with the total number of electrons, they are 
not affected by exp(i?<;), and are, thus, insensitive to the value of 77. We then choose 77 = 1 in order to eliminate the 
interaction of electrons with symmetrical phonons from the transformed Hamiltonian (see Eq. ^ ). We can not do 
the same with 6, since Rb does affect the purely electronic terms. Our approach consists in treating 5 as a variational 
parameter, to be determined by minimizing the total energy. 

The i?h-transformed Fermi operators read 

e^^ci^e-^" = ci, [Ch(7oB) - Sh(7oi3)] , 

e«^C2,e-^-' = C2, [Ch(7oS) + Sh(7oi3)] , (22) 

where we have defined the anti-Hermitian operator B = S{b'^ — b), and Sh(a;), Ch(a;) are the hyperbolic sine and 
cosine functions. Notice the asymmetry with respect to site exchange due to the fact that the antisymmetric phonon 
operator has been defined as a difference of local operators in a well defined order (see Eq. p^ . 

Using Eqs. (p^), we can easily work out the transformation of Hd under Rb- Notice that electron number operators 
remain unchanged, and the only non-invariant terms of Hd are 

(cLci, + i7.c.)e-«- = Ch(27oS) {clc,, + H.c.) + . . . (23) 

and 

e«''(cl^cl^C2xC2T + i?.c.)e-^^^ = Ch{A-ioB){c\^c\^C2^c^^ + H.c.) + . . . , (24) 

where the dots stand for terms that contain Sh(const. x _B), which average to zero in the squeezed-phonon state, as 
discussed below. 

To eliminate the phonon operators, we take the average of the transformed Hamiltonian over the squeezed-phonon 
wave function 

|*sq)=e-"('"'-'''^')|0)ph. (25) 

It is then straightforward to show that 

(^',q|Sh(n7oS)|^',q) - , 



(«'sq|Ch(n7oS)|^',q) = exp 



(26) 



The result is a polaronic effective Hamiltonian iJpoi = Hc\{el,t* ,X* , U* ,V* , J*, J*, P*) + i?ph, where 



el^e-hVL-il[l + 6{2- 5)] , t* = rt , X* ^ tX , U* = U - 2hnj^[l + S{2 - 5)] , 

V* ^ V - 2hn"f^[l - 6{2 - 6)] , P*^T^P, Epi,^hQ[Sh^{2a) + l] , (27) 
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with T = exp [— 25^7Qe~''"] . The parameters Jz and Jxy remain unchanged (thus, they are no longer equal to P*). 
Notice that the energy of the phonons in the ground-state is given by the last term in Eq. ^ Leaving aside hn/2, 
coming from the s-phonons, the contribution of the 6-phonons to this energy can be written as 

El^^lnnCHAa). (28) 

Thus, due to the squeezing effect, we find an increase of the zero-point 6-phonon energy, as if the frequency of the 
latter were J7* = riCh(4Q;). This renormalization effect is not directly due to the electronic interactions, since it does 
not depend explicitly on them. There is, however, aniadirect influence due to a being determined by minimization 
of the total energy, which includes the electronic partO 

The polaronic Hamiltonian is easily diagonalized, and its energy eigenvalues and eigenvectors are listed Jji Table |. 
Leaving aside differences in notation, these results are consistent with those obtained by other authorsEj'tj. 



IV. CORRELATION FUNCTIONS 



The quantities accessible to measurement, besides the energy eigenvalues, are expressed in the form of various 
correlation functions (CF), which can be evaluated through the eigenstates of the effective Hamiltonian, listed in 
Table |. We will indicate by the CF for a generic operator X in one of the eigenstates \i) of the effective 

Hamiltonian. The double average has to be understood as follows: 

^ (^KvI/^qle^'-e^^Xe-^^e-^-l^sq)!^) ■ (29) 

In particular, — E^. We will evaluate a number of correlation functions which have been considered in the 

literature, whose definitions are briefly clarified below. 

The squeezing parameter is directly accessible to experimental measure through the Debye- Waller factor, defined 
as = {{uf)) — ((ui))^. Indeed, by using the time-dependent representation of the displacement through the 

effective Hamiltonian H* , i.e. Ui{t) = exp{iH*t/h)uiexp{—iH*t/?i) one obtains 

pDW = ^[l + e*"] (30) 



where L = ^Jfi/2MQ, is the characteristic plionon length. 

An important quantity is the ratio of the kinetic energy of the interacting system to that on the non-interacting 
one, defined as 

^kin ^ ((^Ea(«L«2a + <Aoaia)))lo^O /g-^N 

" [(^|iE.(«L«2a + «LaiJK>]7o=o 

For iV = 1, F^™ is just the hopping reduction factor r. For N = 2, ^"|^™ vanishes, except in the state \Sb), for which, 
defining 0* like 9 in Table |[ but using the renormalized interactions, we obtain 



^ki„ _ Tsin(2r) 
sin(26l) 



To measure the relation between the charge and the phonon number on each site (related to the polarizability of 
the medium) Ref. ^ introduced, for the case of one electron, the "electron-phonon CF" . In our case it has to be 
defined as 

F;Pee (j,fc = l,2). (33) 

As k can coincide or not with j, we will consider both on-site (j = k) and inter-site (j ^ k) electron-phonon CF's. 
The inter-site charge CF is defined as j{{nin2)), and the longitudinal and transvsrse spin CF's are ((S^f 5*1)) and 
5^)), respectively. The bipolaron "mobility", or singlet- superconductor CF,c2l) describes the hopping of two 
particles occupying the same site, namely, 

{{{c,^c,^c\^cl))) = (vl/,q|Ch[47o(5(6t - 6)]|vI/,q)(^|c2^C2^4^cl J£) 

= r'{l\c,^C2AAl\^) ■ (34) 
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Finally, the (phonon-induced) average charge transfer is defined, following Ref. |T^, as 

^CT ^ (((^i -n2)(&t + &)))„) ^ -26jo{£\{ni - n2)^\£) 

The analytical expressions of various CF's are collected in Table ||. 



V. RESULTS 



The ground state of iJpoi for several values of 50 was identified by searching the minimum of the total, i.e. electronic 
plus phononic, energy upon independent and simultaneous optimization of the parameters that define the shape of 
the orbitals (F), the displaced-oscillator strength (6), and the squeezing-effect strength (a). Given a and S, then 
the optimal hopping reduction factor r = exp(— 2(5^7Qe~^") follows. Specifically, the optimization proceeded in three 
steps. First, the total energy for each eigenvalue was separately optimized with respect to F, S and a. Second, the 
minimal total energy was selected with the corresponding values of the variational parameters. Third, the higher-lying 
eigenvalues where recalculated with F, S and a set by the optimization of the lowest-lying one. This procedure gave 
us the true full spectrum of eigenvalues for each ground state. 

We will now discuss the results for iV = 1 (quarter filled case) and N = 2 (half filled case). 



A. Quarter filled case 



The interplay between the electronic and phononic subsystems is evidenced in Fig. |i| which shows aF, S and exp(— 4a) 
as functions of a, for go = 0.447 eV and hfl = 0.1 eV. The displacement parameter 6 grows smoothly with a, until, at 
a critical length Uc ~ 2.lA, it jumps discontinuously up to (5 ~ 1. The squeezing-related quantity exp(— 4a) decreases 
with a up to Oc, indicating that the squeezing effect becomes stronger to counteract the growing displacement. At 
a = Oc, exp(— 4a) recovers abruptly a rather large value. The orbital-shape parameter aF rises almost linearly with 
a, except around Oc where it suffers a small, sharp increase. Therefore, for = 1, the changes of the phononic 
parameters have only a small effect on F. 

The changes in the electron and phonon subsystems depicted in Fig. |^ are mirrored by the correlation functions 
shown in Fig. ||. We notice that the discontinuities in S and exp(— 4a) at Oc correspond to the localization of 
the polaron, as signalled by the fact that for a = Uc the kinetic energy vanishes (discontinuously) and the on-site 
electron-phonon correlation jumps to a very high value. The localization is driven by r = exp [— 2(5^7oe~'*"] , which 

renormalizes the effective hopping (see Eq.^2^)). The sudden changes observed in Fig.l can be attributed to the 
competition between the squeezing effect, that tends to keep r finite through the factor e~*", and the polaron effect, 
that tends to reduce r through the factor S"^ in the exponent of r. As S increases with a (Fig. |l]), a also increases to 
gain electronic energy by keeping t at a non- vanishing value (see Tabic |) . The price to pay is an increase of the phonon 
energy, Eq.(^). The behavior shown for a < Oc by all quantities in Fig. 2 clearly reflects the continuous variations 
of the phononic parameters of Fig.l. When the phonon energy becomes too large, at a = Oc, the system jumps to a 
state with a small a and zero hopping. However, if the bare hopping amplitude is by itself small (what happens for 
large values of a), the loss in electronic energy when t is reduced is not very important, a remains small, and there is 
no sudden change of state. This can be seen in Fig. ||, where the kinetic energy and the electron-phonon correlation 
functions are plotted as functions of go for two values of a, the large-a case showing a continuous transition. 

The localization critical length depends on the strength of the electron-phonon interaction. We can then draw 
a "phase-diagram" for the system, as shown in Fig. ^. In both panels, the region to the left (right) of each line 
corresponds to itinerant (localized) states. On the left panel we have chosen to plot the bare hopping parameter t, 
calculated at Qc, as a function of go. We also show, on the right, that 5o/^^ ^ good scaling variable for not too 
large dimer lengths, as the curves for different phonon frequencies coincide. For large a, the curves tend to separate, 
enhancing the curvature so as to avoid crossing the x-axis. This implies that the polaronic transition as a increases 
becomes a continuous crossover for small values of go- 



B. Half-filled case 



The half-flUed case presents a more complex behavior than the quarter-flUed one, due to the electron-electron 
interactions that become effective in this case. We first discuss the "phase diagram" , Fig. in order to give an overall 
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view of the physics. Again there are two dominant regimes, locahzed and non-localized, respectively above and below 
the continuous line in Fig. ||. The physics of the system is dominated by the singlet bonding state \Sb) (see Table |), 
which is always the lowest-lying energy level. However, for large a \Sb) becomes degenerate with either the triplet 
states |T) (Fig. ||, top panel) or the charge transfer states \CT) (Fig. |[ middle and bottom panels), depending on the 
value of f/o- These degeneracies just mean that these states are no longer appropriate to describe the system, which 
has become a pair of isolated atoms. For small 5q, as the dimer lenght increases, the system evolves continuously 
from an extended singlet state to a localized state with one electron per site, that occurs for very large values of a. 
For large go a bipolaron is formed (in the \CT) state) which eventually localizes in one of the sites, leaving the other 
one empty. 

While one can draw an fJ-independent phase diagram, the correlation functions do depend on fl. We have therefore 
selected ?ir2 = 0.1 eV, which might be realistic in HTSC's and colossal magnetoresistance materials. We will present 
the results for 5o/^^ = 2.0, 2.2 and 2.5 eV, corresponding to go = 0.447, 0.470 and 0.500 eV, respectively. The 
dotted horizontal lines in Fig. |^ visualize the sections of the phase diagram that we are going to discuss in detail in 
Fig. ^, where the full eigenvalue spectrum is shown as a function of a for the three mentioned values of go- When 
go/^^ = 2.0 eV, 15*6) is always the ground state, with the \T) state approaching it asymptotically for large a (top 
panel of Fig. ||). For g^/hfl — 2.2 eV we find a re-entrant behavior in a range of values of a in which a rearrangement 
of the full spectrum occurs. Finally, in the case 5o/^^ = 2.5 eV one crosses the region of the phase diagram in Fig. 
H bounded by the broken line, where, as we will see below, both the (antiferro) magnetic CF and the charge transfer 
have large values in the 15*6) state. Correspondingly, one sees from the bottom panel of Fig. ^ that inside this CT-AF 
region the state |CT) becomes lower in energy than \T). 

The changes in the electronic and phononic parameters underlying the level crossing in Fig. ^ are shown in Fig. |^. 
In the top panel oF shows an approximately linear trend, steeper for 5o/^^ = 2.0 eV than for go/hfl = 2.5 eV. This 
can be explained by anticipating that, in the latter case, even when, for a < I.60A, the stable state is \Sb), still there 
is an appreciable charge transfer, indicating that the electrons tend to keep apart. Their charge-density clouds can, 
therefore, spread out more largely around the ions (i.e. F can decrease) without too much increase in the Coulomb 
energy. In the intermediate case go/hfl = 2.2 eV, aT jumps between the two lines. The jumps mark the re-entrant 
transition between the localized and the non-localized regime, which therefore implies marked modifications in the 
Wannier functions shape. The middle panel shows the corresponding behavior of the phonon displacement parameter 
6. We see that the passage from the non-localized regime to the localized one is accompanied, as expected, by an 
abrupt jump from (5 ~ to (5 ~ 1. Only for go/fiQ, = 2.5 eV the limiting value ^ = 1 is reached by two jumps 
connected by a range of gradual increase, associated with the crossing of the CT-AF region in Fig. The behavior 
of the squeezing-related parameter exp(— 4a), shown in the bottom panel, is also characterized by sudden transitions 
at the line boundaries. Its dependence on a, combined with the negligible value that 5 takes when 15*6) is the non- 
degenerate ground state, implies the existence of a threshold value for go/^^ below which the electrons influence the 
phonons, but the converse is not true, in the sense that the latter are completely screened out, even though go is not 
small. Large 5 values in the \Sh) state can only be found in the CT-AF region, which, however, is located above the 
threshold. 

In Fig. ^ we show the three largest effective electronic interactions i*, U* and V* (top, middle and bottom panel, 
respectively), for the chosen values of go- Besides showing sharp discontinuities at the critical values of a, their behavior 
clarifies the nature of the CT-AF region of the phase diagram. Indeed, we can see that this region is characterized 
by a significant reduction of the on-site effective Coulomb repulsion and an increase of the inter-site one, which are 
necessary conditions for the stabilization of a local bipolaron, while the effective hopping is only slightly reduced. 
Concomitantly, the bipolaron mobility remains finite (Fig. ^, bottom panel) while charge-transfer fluctuations build 
up (Fig. [1^, top panel), and the inter-site spin (or magnetic) correlations are still important (Fig. bottom panel). 
We would like to remark that the existence of this region of the phase diagram (probably the most interesting one 
as far as high- Tc superconductors are concerned) depends strongly on the squeezing of the phonon states, as can 
be seen in the bottom panel of Fig. 0. In other words, the squeezing effect counteracts the oscillator displacement, 
and succeeds in preserving a rather itinerant character even in the presence of both appreciable charge transfer and 
magnetic correlations. 

VI. CONCLUSIONS 

We have presented here a variational study of a dimer where interacting electrons occupy a single orbital per site and 
are coupled to Holstein phonons. We have put into evidence that the variation of the dimer length and/or the strength 
of the electron-phonon coupling yield a strong renormalization of the interaction constants as well as the width of the 
Wannier functions describing the local orbitals, which establishes a link between the electronic interactions and the 
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phonons at a deeper level than predicted by the standard polaron approach. The squeezing of the phonon states has 
been shown to be particularly relevant, yielding sharp changes in the state of the system as the dimer length is varied 
at fixed electron-phonon coupling. The experimental interest of this finding is that situations where the interatomic 
distances can be varied, even by a small amount, either uniformly by an external pressure, or inhomogeneously by 
dimerization in a chain compound, can lead to important changes in the system's behavior. 

Obviously one has to be careful in making predictions about the behavior of an extended system based on calcu- 
lations for a two-site cluster. Nevertheless, despite the fact that Figs. ^ and |^ cannot be viewed as phase diagrams, 
we believe that they allow us to speculate on the nature of the corresponding states in a lattice. Remaining in 
the half-filled case, the region where the extended bonding singlet is stable should show the antiferromagnetic order 
characteristic of the usual Hubbard model, with a continuous evolution towards a local-moment paramagnet as the 
interatomic distance is increased. For a doped system, one could expect the two-site singlets to build up a RVB-like 
state, or an inter-site bipolaron, in the lattice. On the other hand, the |CT) state should correspond to a charge 
density wave in the lattice. In connection to the HTSC materials, our most significant finding is that there is a region 
in the parameter space where the kinetic energy, the charge transfer and the magnetic correlations are simultaneously 
large, while the phonon state is partially displaced and strongly squeezed. This result suggests that it is perhaps not 
necessary to set a sharp alternative between either magnetic or charge instabilities. Based on our results, one might 
speculate that, in the normal state, those materials might be in a situation corresponding to the one in the dimer, 
where both types of fluctuations can be simultaneously large. It is interesting to notice that the existence of this 
narrow region (between the dashed and continuous lines in Fig. ^ in our model is entirely due to the squeezing of 
the phonon states. Indeed, we checked that it disappears if one keeps the squeezing-parameter a fixed at zero, giving 
rise to a direct transition from the extended singlet ground-state {S « 0.) to the charge transfer one {6 = 1.). We 
thus suggest that the competition between magnetic and charge instabilities in HTSC materials might be driven by 
electron-phonon interactions in the presence of strongly squeezed phonon states. It would be interesting to use the 
dimer solutions obtained here as building blocks for lattice states, with the aim of obtaining quantitative support for 
our qualitative analysis of the lattice case. Work on this line is now in progress. 
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TABLE I. Eigenvalues and eigenvectors of the Hamiltonian of Eq.(l) for A'^ = 1, 2, 3, 4. D is the degeneracy of the state. The 
labels a, b indicate bonding and antibonding character, while S and T correspond to Sing let and Triplet states. For N = 2 we 
have dcfinod Eu = 2eo + U + P, Ev = 2eo + V + J^y, r = ^{Eu - EvY + 16(t - Xy, tan(^) = -A{t - X)/(Eu - Ev + r) 
and T = t - 2A. 



Filling and Energy 


D 


s, 


Eigenvectors 


(Ar = 1) 

£2 = eo + f 


2 
2 


1/2, 1/2 
1/2,-1/2 
1/2,1/2 
1/2,-1/2 


|16,T) = ^[ct^-Kc+^]|0> 
11a, T> = -i=[ct^ -ct^]|0) 


(AT = 2) 

Esb = ^{Eu + Ev -r) 


1 


0,0 




Et,±i = 2eo + V - ,7^ 
Et,o = 2to + V - J^y 


2 
1 


1>1(-1) 
1,0 




EcT = 2u) + U ~ P 


1 


0,0 


\CT) = ^[cici^ - ci4^m 


Esa = iiEu + Ev+r) 


1 


0,0 


\Sa) = ^[cose(ct^ct^ -^4^4^) 



(JV = 3) 

Ssb, = 3eo + U + 2V-J,-T 
E3a =3eo + U + 2V - J^+T 



1/2,1/2 
1/2,-1/2 
1/2,1/2 
1/2,-1/2 



136, T) = ^14^014^+ 4^4^c\^]\0) 

\3bA) = :}^[4,4^4^ + 4,4^4^]\o) 

|3a,l> = ^[ct^ct^ct,-ct^ct^ctj|0> 



(JV = 4) 








Ei = 4eo +2U + 4V - 2J^ 


1 


0,0 


|4> = ct^cl^ct^ctjO) 
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TABLE II. Analytic expressions of some correlation functions ((-X')) in the indicated eigenstates. 



X 


IM 


\Sb} 


|T,±1) 


|T,0> 


\CT) 




1 

4 


sinh^(2a) 




isinh2(2a) + 7„2 


i sinh^(2a) 


isinh2(2a) 


i sinh2(2a) 


u — ^> 


+7?(l + «)' 






+To 


+7o 


+7o U + 


^jl>lb^ 


1 

4 


sinh^(2a) 




isinh^(2a) +7^ 


i sinh^(2a) 


isinh^(2a) 


i sinh^(2a) 


{jjik = 1, 2) 


+7?(1 - Sf 




-7oi5(2 - (5) sin'' 6 


+7o 


+7o 


+7o(l-'5)^ 


1"1"2 


0. 


i cos^ e 


1 
4 


1 
4 


0. 


Si S2 


0. 


-icos^e 


1 
4 


-icos^e 


0. 




0. 


-|cos^e 


0. 


1 

2 


0. 




0. 




0. 


0. 


2^ 


(m -n2)(6t -6) 


-2570 


— 16(570 sin^ $ 


0. 


0. 


— 16570 


^/mf-br)) 


^cosh(4a) + (2570 )2 


cosh(4a) + (4i57o sin Oy- 






^cosh(4a) + (4570)2 



10 



FIG. 1. The Wannier-function-shape parameter aT, the displaced-oscillator parameter 5, and the squeezing parameter 
exp(— 4q) as functions of a, for A'^ = 1, hQ, — 0.1 eV, and go = 0.447 eV. 



FIG. 2. Variation of the reduced kinetic energy, the charge transfer, and the on-site and inter-site EP CF's with a, for A'' = 1, 
no, = O.leV and go = 0.447 eV. 



FIG. 3. The same quantities as in Fig. g as functions of go, for fixed length values a = 2.3 A and a = 4.0 A. The inter-site 
EP CF for a = 4.0 A has vanishingly small values on the scale of the figure. 

FIG. 4. The "phase diagram" for N — 1 and two values of hi^. The vertical coordinate is the bare hopping parameter t, 
and the horizontal one is the electron- phonon coupling go (left panel). The electron is itinerant on the left of each curve, and 
localized on the right. The right panel shows that go/hQ is a good scaling variable for small a, as the two curves coincide. 

FIG. 5. The "phase diagram" for N — 2. The electrons are non- localized below the full line, and localized in the region above 
it. The broken line delimits the region where both the magnetic CF and the CT are large. The horizontal lines correspond (for 
hii = O.leV) to the values of go/hfl considered in the numerical analysis. 



FIG. 6. Dependence of the full eigenvalue spectrum for N = 2 with a, for fifi = 0.1 eV, and the indicated values of go- 

FIG. 7. Microscopic electronic and phononic parameters for A'^ = 2 vs. a for the indicated values of go- Top: the Wannier 
wavefunction shape parameter aP; middle: the oscillator displacement 5; bottom: the squeezing parameter exp(— 4q). 

FIG. 8. The effective hopping parameter t* (top), the on-site interaction U* (middle), and the inter-site interaction V* 
(bottom) as functions of a, for N = 2, TiQ. = 0.1 eV, and the indicated values of go- 



FIG. 9. The on-site electron-phonon CF (top), the inter-site electron-phonon CF (middle), and the singlet-superconductor 
CF (bottom) as functions of a, for N — 2, h^l = 0.1 eV the indicated values of go- 



FIG. 10. The charge transfer (top) and the absolute value of the magnetic CF, |((S'i5'2))|, (bottom) vs. a, for N = 2, 
hO, — 0-1 eV, and the indicated values of go- 
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